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Abstract 

It is well known that a spherically symmetric constant density static star, 
modeled as a perfect fluid, possesses a bound on its mass m by its radial size 
R given by 2m /R < 8/9 and that this bound continues to hold when the 
energy density decreases monotonically. The existence of such a bound is in- 
triguing because it occurs well before the appearance of an apparent horizon 
at m = R/2. However, the assumptions made are extremely restrictive. They 
do not hold in a simple soap bubble and they certainly do not approximate 
any known topologically stable field configuration. In addition, such configu- 
rations will not generally be compact. We show that the 8/9 bound is robust 
by relaxing these assumptions. If the density is monotonically decreasing and 
the tangential stress is less than or equal to the radial stress we show that 
the 8/9 bound continues to hold through the entire bulk if m is replaced by 
the quasi-local mass. If the tangential stress is allowed to exceed the radial 
stress and/or the density is not monotonic we cannot recover the 8/9 bound. 
However, we can show that 2m/ R remains strictly bounded away from unity 
by constructing an explicit upper bound which depends only on the ratio of 
the stresses and the variation of the density. 
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I. INTRODUCTION 



Consider any static solution of the Einstein equations with the matter satisfying the null 
energy condition. The Penrose singularity theorem shows that this system cannot have a 
trapped surface. In a spherically symmetric configuration, the first apparent horizon in the 
initial data occurs when the ratio of (twice) the mass to its radial extent, 2m/ R, is one. 
However, it is also well known that for a spherical star composed of ordinary matter with 
positive energy described as a perfect fluid with a monotonically decreasing energy profile, 
2m/R cannot exceed 8/9, the constant density value 0-0]. Such a bound is particularly 
interesting because it occurs strictly before the appearance of an apparent horizon. Here we 
wish to examine this bound more closely when the assumptions underlying it are relaxed. 
Several very strong assumptions on the distribution of matter enter its derivation. Even 
in an astrophysical stellar object where matter is described phenomenologically, it is not 
clear that the perfect fluid assumption is justified. A humble soap bubble consisting of 
a membrane with a given (tangential) surface tension supported by the pressure of the 
enclosed perfect gas violates both the monotonicity of the energy density and the perfect 
fluid assumption. The approximation clearly also does not represent accurately the interior 
of topological defects such as monopoles || . The balance of forces providing the equilibrium 
typically turns out to be analogous to that holding a soap bubble together. In addition, 
in many extensions of Einstein gravity, the effective stress tensor describing a perfect fluid 
does not assume the perfect fluid form. Field configurations typically will not be compact, 
in which situation we require a generalization of the mass that holds thoughout the bulk. 
This requires the replacement of m by a quasi-local mass. 

While the 8/9 bound is not a universal one, it is robust in the sense that under physical 
conditions which are at least reasonable classically, the mass continues to be bounded by a 
value strictly below the apparent horizon value, m = R/2. It appears to be impossible, even 
in principle, to construct a static distribution which saturates it. 

In Sections II and III, we establish our notation. We show in Section III that if the 
matter satisfies p + Sc > 0, where p is the energy density and Sc is the radial stress, 
the object cannot have an apparent horizon and thus 2m/ R is strictly bounded away from 
1. This condition, p + Sc > is one of the so-called 'null energy' or 'null convergence' 
conditions |7|-[J • It is interesting because we need no restriction on the tangential stress nor 
do we need to assume p > 0. In a spherically symmetric geometry, the tangential stress 
will generally differ from the radial one except at the center where the constraints dictate 
that they coincide. Consider the ratio, 7, of tangential to radial stress. In a perfect fluid 
7 = 1. We further show, again in Section III, that static matter satisfying p > together 
with 7 < 1 must have positive radial pressure which monotonically decreases outwards. This 
guarantees that p + Sc > which provides another way of excluding apparent horizons. 

The rest of the article is devoted to investigating how close 2m/ R can get to one, and 
generalizing the 2m/ R < 8/9 bound noted above. We summarize very briefly the simple 
constant density star in Section IV. In Section V we consider 'stars' that are monotonic with 
positive density. If 7 < 1, it is simple to show that the 8/9 inequality continues to hold, not 
only on the boundary but through the entire bulk. Indeed, the configuration need not be 
compact. If 7 > 1 anywhere, however, we obtain a slightly weaker result. We construct a 
bound which shows that 2m/ R is strictly bounded away from unity. If 7 TOaa; approaches 1 the 
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bound smoothly approaches 8/9; as 7 becomes unboundedly large the bound approaches 1. 
In particular, for a monotonic star with positive radial pressure and for which the transverse 
pressure is less than the density we can show that 2m/ R < 0.974. 

Finally, in Section VI, we relax the assumption of monotonicity and find essentially the 
same results, except that now the bound depends both on the variation of the matter as 
well as on j max . 



II. STATIC LIMIT OF EINSTEIN EQUATIONS 

The spacetime metric describing a static solution of the Einstein equations can always 
be written in the form 

ds 2 = -N 2 dt 2 + g ab dx a dx b , (1) 

where N is the lapse function and the shift vanishes. N is also the norm of the global timelike 
killing vector, dt, and so must satisfy N > 0. The spatial geometry at constant t is described 
by the metric tensor g ab . Both the material current vector J a and the extrinsic curvature 
tensor K ab (describing the embedding of a hypersurface of fixed t in spacetime) vanish. In 
the canonical formulation of the theory, the momentum constraints of the theory are then 
vacuous. The hamiltonian constraint reduces to the form |L0|,nH (see also the appendix to 



■R=16np, (2) 

where 71 is the scalar curvature constructed with the spatial metric g ab and p is the material 
energy density. Given some specification of p, Eq.([|) is a constraint on the spatial geometry, 
g a b. It does not involve the stresses operating on p. In the spherically symmetric case we 
will see that the intrinsic geometry is completely specified by p. The advantage of working 
within the canonical formulation is that this constraint is isolated explicitly. 

Given that the time direction is Killing the evolution in this direction must be trivial. 
The dynamical Einstein equation reduces to K ab = 0, and now reads 

- V a V b N + 7Z ab N = 8tcN (S ab - -g ab ti S + -<? o6 p) , (3) 

where V a is the covariant derivative compatible with g ab , lZ ab is the associated Ricci tensor, 
S ab is the material pressure tensor and S is its trace. In a perfect fluid the stress is isotropic 
with S ab = Pg ab - The other evolution equation, g ab = 0, is trivially satisfied. 

For given p and g ab , Eqs.([|) consist of six PDEs for the seven functions, N and S ab . This 
counting is not very precise, because if we had a realistic fluid/field theoretical model we 
would have to supplement these equations with an 'equation of state' which would convert 
the equations from an underdetermined to an overdetermined system. 

We will suppose that the spatial topology is R 3 . For an object with energy density of 
compact support (a star) or falling off sufficiently rapidly at infinity the spacetime will be 
asymptotically flat with iV — > 1 at infinity. The appropriate boundary condition on S ab in 
an object of compact support is that its normal component vanishes on the surface. For 
much of our discussion these boundary conditions are irrelevant. 
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Taking the trace of the equations, (|3|), and eliminating 1Z in favor of p using Eq.(@), we 
obtain the linear elliptic equation for N, 



AN = 4n(p + tr S)N . (4) 

If the strong energy condition is satisfied we have p + tr S > everywhere, and so AN > 
when N > 0. Thus the solution cannot possess an interior maximum. N falls towards the 
center. Even if the potential, V = p + trS, is large, N never fall to zero in the interior (e.g., 
|T2|1). The lapse can go to zero only when the density or pressure becomes unboundedly 
large. 

Let us assume that N is positive in the exterior and negative on a compact region W. 
N vanishes on dW but will have positive outward gradient. If we integrate AiV over W we 
can turn it into a surface integral which must be positive. On the other hand, from Eq.([|) 
we see that AiV < on W, so we have a contradiction. We also see that it is impossible for 
N to just touch zero at a point. At that point we would have that N, its first derivatives 
and its second derivatives all vanish. Thus the function could never grow away from zero 

One can deduce the conservation law, 

V b S ab = -(S ab + pg ab )^f } (5) 

directly from the static Einstein equations, Eqs.(^) and ([!]). To do this, we simply take the 
divergence of Eq.(|3|). Exploiting the Ricci identities, 

[V a ,V b ]V b = R ab V b , (6) 

and the contracted Bianchi identity for lZ ab , V a TZ ab = V b 7Z/2, we reproduce Eq.(||). 

It is clear from Eq.(ffl) that there are no non-trivial vacuum static solutions in the theory. 
We have AiV = everywhere. If there is no internal boundary, the solution is iV = 1 
everywhere. Now lZ ab = 0, as well as TZ = 0, so that the geometry is flat everywhere. There 
is a well known result that the only perfect fluid static equilibria are spherically symmetric 
|14] ) [15| . This result implies that for perfect fluids the spherically symmetric analysis is 
complete. A discussion of the symmetries of equilibrium configurations is provided in |[16| . 



III. SPHERICAL SYMMETRY 

The line element describing the spatial part of a spherically symmetric geometry can 
always be written as 

ds 2 = d£ 2 + R 2 dQ 2 , (7) 

£ is the proper radial distance on the hypersurface, R is the areal radius. For R 3 topology, 
£ has domain [0, oo). The appropriate boundary conditions on R are 

R(0) = 0, dR/d£\ = R'(0) = 1 . (8) 

The scalar curvature TZ is given by 
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2 {RE!) - R 12 -I , (9) 



where primes denote derivatives with respect to £. The constraint equation can be cast in 
the form (see, for example, |T7| , P^] ) 

R' 2 = i- 2 -^, (io) 

where the positive quasi-local mass is given by 

m = 4vr [ £ pR 2 R'd£ = 4vr [* P R 2 dR. (11) 
Jo Jo 

It is immediately clear that 

m<R/2 (12) 

everywhere. In general, R' 2 < 1 in any regular geometry when the weak energy condition 
(p > 0) is satisfied, so that R < £ everywhere flI8|. To show this we substitute Eq.(pf) into 



Eq.(|) to get 

2RR" + R' 2 - 1 = -8nR 2 p . (13) 

At the center we have R' — 1 and m — 0. From Eq.(|ll]) we see that m increases as soon 
as we meet matter and thus R' drops below 1. Let us assume that it later rises up to +1. 
However, from Eq.flTB]) we see that at this point R" < so it cannot be rising! On the other 
hand R' can drop below —1. We again get R" < which means that it cannot ever rise 
up again to the asymptotic R' ~ 1. Thus for any regular spherical geometry satisfying the 
weak energy condition — 1 < R' < +1 and R' = +1 only at the origin and at infinity This 
holds true for any solution of Eq. (|T3"D , no static assumption is required. 



It is clear from Eq. ( JTOj) that m is positive everywhere in a regular geometry and vanishes 
only at the center and in any vacuum region surrounding it. In a static configuration, the 
extrinsic curvature vanishes so that an apparent horizon is a minimal surface with R' = 0. 
Thus, if the geometry is free of an apparent horizon, it must also be free of singularities. In 
such a geometry < R' < 1 and m increases monotonically with I (or R). We emphasise 
that the spatial geometry and with it the ADM mass is completely determined by the source 
energy density. The material stresses play no role whatever. 

At the surface of a compact object of radius R = Rq, the quasi- lcoal mass coincides with 
the constant ADM mass, mo- The exterior solution is given by Eq.flTOD 

*- = l-^. (14) 

In a spherically symmetric geometry, any symmetric tensor is completely characterized 
by two scalar s. We have 

Tlab = R-cn a nb + R-R^gab - n a n b ) 

Sab = S c n a n b + S R (g ab - n a n b ) . (15) 
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Here n a is the outward pointing normal to a two sphere of fixed proper radius. The two 
scalars appearing in the Ricci tensor can be expressed in terms of 7Z, R and R' as follows 
0: ' 

Kc = \n - j- 2 (i - R") 

Kr = \k + ±- 2 (1-R' 2 ). (16) 

Taking the two independent projections of Eq.(^), we therefore have in any spherically 
symmetric static equilibrium, 

N" = N {47r(p - S c + 2S R ) - ^} (17) 
R'N' = RNU-kS c + ^\ . (18) 



We can also combine Eqs.flT8|) and (]13[) to obtain 

Till p/ AT> 

-B + BN = *^ + 8 ^ (19) 

These three equations are the complete set of equations satisfied by any static spherically 
symmetric system. 

If the matter satisfies 



P + S c > , (20) 

we can immediately deduce from Eq.([l9|) that any spherical static configuration cannot 
have an apparent horizon. The apparent horizon coincides with a minimal surface, i.e., 
when R' = 0. From ( |T9"D at a minimal surface we see that R" < 0. However, at the 
outermost minimal surface we must have R" > since the area is increasing outwards. 
Thus we have a contradiction. Thus we have shown that, in a static star satisfying Eq. (p0|) , 
R' > 0. From Eq . (|10"D this is equivalent to showing m < R/2. In || essentially the converse 
of this argument was given, showing that if a minimal surface existed (the throat of a static 
wormhole) then the matter cannot satisfy Eq.(|20D. This result has been recently proven in 
|i9| . However, that proof requires that both p and Sc be positive while we only need to 
impose a condition on the combination. 

The energy condition Eq. (pOf ) is a single component of what is called the 'null conver- 
gence' or 'null energy' condition 0,0]. It is a consequence of each of the three standard 
energy conditions (the 'strong', 'weak', and 'dominant'). Consider the outgoing radial null 
vector 

e = (l/N, 1,0,0), (21) 
and multiply it into the spacetime Ricci tensor ^TZ^ to get 

(4) 7^er = {A) ever = Mp + s c ) . (22) 
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where ^Gfw is the spacetime Einstein tensor. The equality above is to be expected because 
the Einstein tensor only differs from the Ricci tensor by a trace and the trace term, when 
dotted twice with a null vector, vanishes. The positivity of ^IZ^^F, 1 ' implies Eg . (p0|) . 
Choosing £ M to be an outgoing tangential null vector we obtain p + Sc + 2Sr > 0. 



If both Sc > and p > hold independently, as supposed in it is clear that Eq. pOD 
is satisfied which guarantees R' > 0. We also now have that the right hand side of Eq.(|l^) is 
positive so that N' > everywhere. The lapse function, the length of the Killing vector, for 
any regular solution must grow monotonically out from the centre to its asymptotic value 
one. With positive radial stress and positive p we do not need to assume the strong energy 
condition. Spherical symmetry is very restrictive. Compare this to the spherically symmetric 
statement of the maximum principle which was applied earlier to the trace equation, Eq.(|j). 

For completeness, we note that the lapse is evaluated in the exterior of a compact object 
as follows: using Eq. (|IBD , we have 

R ' N > = N ^ , (23) 
R A 



£ = -2=0". (24) 



so that using Eq.fll 



R' R 2 \ R 
The boundary condition at infinity, N — > 1, fixes 

Eq.(p5D together with Eq. flTj]) reproduce the exterior Schwarzschild form of the spacetime 
metric. 

The conservation of the stress tensor reduces to the single equation, 

R' N' 
S'c + 2—{S £ - S R ) = -(Sc + p)— . (26) 

We deduce immediately that at £ = in a non-singular geometry 

Sc = S R . (27) 

The perfect fluid form of the stress tensor is the only one consistent with the symmetry at 
the origin. This is exactly as in newtonian theory. 

While we only needed a condition on the radial stress to eliminate apparent horizons, 
the transverse stress does play a role in establishing the equilibrium. This can be seen 
in simple mechanical models. For example, in a soap bubble, the surface tension, which 
is effectively a negative transverse stress, is the object which balances the positive outward 
pressure difference between the inside and outside. On the other hand, if we had an evacuated 
spherical metal shell, with a vacuum inside and positive pressure outside, the outside pressure 
forces the metal shell to contract setting up a positive transverse stress. The radial stress 
obviously increases outwards and is balanced by the positive transverse stress. In a self- 
gravitating system, we expect the radial pressure to decrease outwards. However, it need 
not if there are large positive transverse stresses to support the external pressure. 
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In a spherically symmetric geometry it is possible to exploit the first order Einstein 
equation, Eq. flTlf ), to reduce the dependence on the stress tensor appearing in Eq (|l7D to a 
dependence on the ratio of the tangential to the radial stress, 

Sr 

We have 

N" + E (1 - 27) N> = {4vrp - ^ (1 + 2 7 )} N 

In particular, if the stress is isotropic then 7 = 1 and Eq.(^) is independent of S a b. 

Alternatively, we can exploit Eq.(|T8|) to eliminate the lapse from the conservation equa- 
tion, Eq.(|6]) 

S'c + 2^(Sc - S R ) = -£(S C + P) (^S c + ^) . (30) 

In the isotropic limit, Eq.(^) is the Tolman-Oppenheimer-Volkov equation. In the newto- 
nian limit, the RHS of Eq.([30|) is replaced by —pm/r 2 . 

In the simple mechanical models provided above, it is clear that positive transverse 
stresses reduce the internal pressure and vice versa. We can exploit Eq . (|26|) (or Eq.(|30D) to 
show that, in general, if the transverse pressure is smaller than the radial pressure the radial 
pressure builds up inside. We will give two slightly different versions of this result. 

First, let us assume Sc — Sr > 0, p + Sc > and p + Sc + 2Sr > (both of the latter 
coming from the 'strong energy' condition). From the trace equation, Eq.(f|), we have that 
N' > and from Eq.(|l9|) we have R! > 0. When these are substituted into Eq.(^) we get 
S' c < so the pressure monotonically increases inward. The object does not need to be 
compact. 

Alternatively, even more simply, let us assume Sc — Sr > and p > 0. The hamiltonian 
constraint guarantees m > 0. Suppose first that the object is compact. We cannot have an 
apparent horizon outside so we have R' > on the boundary. At the boundary of a compact 
object Sc = and so from Eq.(^) we get S' c < so that Sc is decreasing outwards and 
so must be positive near the boundary. However a positive Sc makes the right hand side of 
Eq.flSOD even more negative so Sc becomes ever larger as one moves inwards. Thus we have 
shown that Sc > and monotonically decreases as one travels out. In turn this guarantees 
both p + Sc > and AtcS c + m/R 3 > 0. Hence R! > and N' > 0. Note that we need not 
assume that Sr vanishes on the boundary but it cannot be positive there. In other words, 
surface tension is good. If the object is nor compact, the argument we have just presented 
is valid in the region bounded by any sphere with Sc > 0. 

In this section, using very weak assumptions, we have demonstrated that in a spherical 
static star we have < R' < 1. From Eq.([3I]) we now get that both m > and 2m/ R < 1. 
However, we have to date no information on how close R' can get to zero, how close 2m /R can 
get to 1. This will be discussed in the next sections, where, by imposing various restrictions 
on the matter, we get extra control on the behaviour of 2m/ R. 



S 



(28) 



(29) 



IV. CONSTANT DENSITY PERFECT FLUID STAR 



It is our good fortune that for a perfect fluid constant density star, Eq.(|3tJ) is exactly 
solvable. Eqs.fllOD and ([11]) reduce to 

R- + (^p) R 2 = 1 . (31) 



We then have 



<jP 47T - g- D2 (P + Po)fp + ^o) , -32; 



dR 1 - f p # 2 

with the well known solution, 



M 2 f* 2m \V 2 



PlR) '* I 7 / 2 tf y M/» I ^ (33) 

The pressure always exceeds the newtonian value. In fact, in an isotropic uniform newtonian 
fluid ball of radius Rq, the central pressure is given by P = P c , where 

P c = y P 2 X ■ (34) 

The pressure given by Eq. (|33"D diverges at the center R = when m = 4i? /9. This occurs 
when the surface lapse, AT = 1/3. If m > 4_R /9, it diverges at some finite value of R. As 
2mo is increased up to Ro, the divergence moves out to Rq. 

Let us now examine the lapse. In a constant density perfect fluid, Eq.fl29|) assumes the 
very simple form, 

(0-o. 

We exploit the continuity of the lapse and its first derivative across Ro which follow from 
Eq.fll7l) and fll8| ) We first integrate out from some interior point to the surface at Ro = R(£o)'- 

(*!.) =(—) =^ (36) 

\ R /R<R \ R /R=R 



where we have exploited Eq.(|25D to evaluate the RHS. We integrate again over the same 
domain. We find for the surface lapse, 

N = N C + ^ f e °d£R(£), (37) 
R'o Jo 

where N c is the value of the lapse at the center. We note that generally 

/Ro / 2m \ -1 / 2 

^ RdR(l-—\ . (38) 



Thus, in a constant density star, 



^.,>( 1 -^.i( 1 -( 1 -!fn5, (39 , 



and we get 



-Ro / 2 \ \ R 

We require N c > 0. Eq.(EO) then implies 



„ 0={1 -^)" = „. + Ul-(l-?Pf). (40, 



m <^Ro, (42) 



or 



exactly as before. This route, however, has the advantage that Eq.(|29|) is linear in N unlike 
Eq.(^0]) which is nonlinear in Sc- 

It is worth noting that iV — > as m — > 4R/9 should not be viewed as the Killing vector 
going null. It is another version of the 'collapse of the lapse' phenomenon, in this case driven 
by the fact that the pressure is becoming unboundedly large. 



V. MONOTONIC STARS 

Buchdahl || demonstrated that if the energy density profile in a star is monotonically 
decreasing, and it is modeled as a perfect fluid, this 4/9 bound continues to hold. The 
constant density star saturates the bound within this class of systems. In this section we 
follow Buchdahl in only considering objects with monotonically decreasing densities but we 
will push the calculations much further. We start off with a perfect fluid assumption and 
rederive the 4/9 bound. We then weaken this to the dominant radial pressure assumption 
($£ > Sr) that we used in Section III and prove that the 4/9 bound is still valid. We next 
extend the inequality to interior points. We finally consider the situation where Sr may 
be larger than Sc- We no longer can recover the 4/9 bound; however if the ratio of the 
pressures is bounded we show that m/Ris strictly bounded away from 1/2. 

Let us define 

4tt . . m . . 

T W = ff . («) 

so that 

is an average of p(R) (not to be confused with the physical average) within a euclidean ball. 
Thus if p' < 0, it is clear that (p) > p and 
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y=T w£o - <45) 



In particular, one can deduce that m/R > uiqR 2 /Rq, so that 

^r 1 / 2 / 2m R 2 \- 1/2 



rj -V-^r) • (46) 

a lower bound is always provided by a constant density star with the same m and Rq. 

We mimic the constant density star calculation. This is essentially the Buchdahl deriva- 
tion, however we allow for a non-perfect fluid. We combine Eqs.(|T7|) and (|l8| ) to give 




N" N'R' 47r N 

- = 1(P ~ (P» + 2 (5« - S £ )] (47) 

Both terms on the RHS of Eq. ([47l) are negative when p' < and Sc — Sr > 0. Thus we 
have 

/ AT>\ ' 

(48) 

with equality only in a constant density star supported by isotropic pressure. The remainder 
of the calculation in this case mimics that for a constant density star. 

As before, we first integrate Eq. (|48|) out from some interior point to the surface at 

R = Rq'. 





(49) 

o JX o 

We follow this by integrating out from the center at i = to the surface. We find 

N >N C + ^ f e °d£R(£). (50) 
Rq Jo 

We require a lower bound on the integral appearing in the second term on the RHS. We 
cast it, as before, in the form (^8|). Using Eq.(f46D, it is clear that 



When we substitute Eq. (^) and Eq. ( [25] ) into ( |5*UD together with the requirement that N c > 
we recover Eq.(|4~l"D and so we have that 2m /i?o < 8/9. 

This gives only a bound at the boundary. If the configuration has a 'thin' atmosphere 
with 3p < (p), m/R is decreasing so the maximum value of m/R occurs somewhere in the 
interior and not on the boundary. In such a scenario the above result is not very useful. 
Happily, the argument can be tweaked to show that 2m/ R < 8/9 through the whole system. 

Let us assume p > 0, p' < and Sc > Sr. The argument at the end of Section III shows 
us that Sc > 0. From Eq.(|l^) we have 
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N' N 



R R! 



m 



— = — 4* Sr. + > — U-r . 



R 3 



N ( m\ 



Rf \R 3 J 



(52) 



Let us assume that 2m/ R possesses a maximum at a point a distance i\ from the center. 
The monotonicity of N'/R (Eq.©) and Eq.fljjD gives (contrast Eq.(H)) 



iV 



1? 



N m\ 



Integrate from the center to l\ to get 

iVi > iV c - 



,iVm\ 



' N m 
> N r + — — 



i2' 



(1 - 2m/i2) 
«i RdR 



1/2 



2m 1 R 2 /Rl) 



1/2 



(53) 

(54) 
(55) 

(56) 



where the last line follows from the monotonicity of m/R 3 . This can be integrated to give 



N!>N C + 



m\Ni 



R\ 



2mA 



1/2 



(57) 



Requiring that iV c > allows us to cancel the N\ on both sides and we immediately get that 
2m/ R < 8/9. 

To deal with the situation where Sc can be less than Sr we need a somewhat more 
complicated argument. We add now as one of our assumptions that Sc > 0. 
If we divide Eq. ( fT7| ) by Eq. (|T8|) we can get 



'W\'_N'R! (p-(p))+2(S R -S £ ) 



R R 



Sc + f 



(58) 



Let us assume that the term on the right hand side of Eq. 
is bounded. In other words we assume 



(p-(p)) + 2(S R -Sc) 
Sc + 4 



which depends on the sources 



(59) 



It is clear that (3 cannot be negative because the numerator vanishes at the center. We have 
(3 = for a monotonic star with Sr < Sc- In general, it will be some positive number. 
Eq.(|58|) now reads 



N' R' 
~R ~R 



(60) 



Find the point where 2m/ R is a maximum (call it £% as before) and integrate Eq.(^) out 
to it to give 
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(61) 



so that 




(62) 



As before, we integrate this equation from the center out to i\ to get 





1 / R 



RJ 1 Jo \Ri 



m 



Rdl 



/ R\P 



o \Rx 



«i ( R 



Rdi 



RdR 



R'R?JiJo \RJ (l-2m/R) 1/2 
N m\ r R i ( R\ fi RdR 



Ft" ft'' J | Jo \nj ( x 2m 1 B^ l/2 • 



(63) 
(64) 
(65) 

(66) 



In going from Eq.(^3|) to Eq.(|64"D we use Eq.flTBD and in going from Eq.(|64"D to Eq.(|65|) we 
use that Sc(£x) > 0. It is clear that the integral in Eq. (|66|) is finite and well behaved for 
any finite (3. Thus we get a bound on 2m /R which is strictly bounded away from 1. Only 
in the limit as (3 — > oo does the integral go to zero. In this case the bound on 2m /R — > 1. 
In the other limit, when j3 — > 0, we recover Eq.(|5*6|) and so we get 2m /R — > 8/9. In the 
special cases where (3 = 2, 4, 6, . . . the integral in Eq.(^) can be done simply. This includes 
one especially interesting case. 

Let us assume we are given a monotonic star with positive radial pressure (these as- 
sumptions can be justified by stability criteria). Let us further assume that the transverse 
pressure is bounded. More precisely let us assume Sr < p. This can be justified on some 
kind of speed of sound argument. From the monotonicity we get Sr < (p). From these we 
immediately get (3 < 6! Now we can do the integration and get 2m/ R < 0.974. 

Alternatively, if the material was approximately a perfect fluid we could use the ratio of 
the pressures, 7, that we introduced earlier. It is clear that (3 < 2{^ Max — 1). 



VI. COMPLETELY GENERAL SPHERICAL CONFIGURATION 

Let us now consider a general static spherical ball. We no longer wish to assume either 
monotonicity or a perfect fluid. The only constraints we place are that both p > and 
Sc > 0. We also assume that (3 as defined by Eq.(^) exists. We are less interested in 
obtaining the tightest bound on 2m/ R than in establishing that such a bound exists. 

All the equations, starting from Eq.(|59| up to and including Eq.(|65"D continue to hold. 
However, in going from Eq.flBSD to Eq.(^Bj) we used the monotonicity. One way of avoiding 
that difficulty is by replacing Eq . (|66D with 
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This uses that 2m/R = 8ir(p)R 2 /3 > 8n(p) Min R 2 /3. Eq.(|67|) can be simplified by introduc- 
ing a new variable x 2 = 87r(p) Min i? 2 /3 and x\ = 87r(p)Min-Ri/3 = 2rh 1 /R 1 < 2m 1 /Ri where 
4vr(p) M in-R 3 /3 < m. We then get from N c > 



m 



> 



(P) 



2rrii 



Ri 2(p) Min Jo 



j-i 



. ' i uO 
•i'l 



(68) 



It is clear that the right hand side of Eq.fl6"8|) is finite and bounded away from zero as long 
as (p)rviin is non-zero. This is very misleading because of the dependence of X\ on (p)Min- If 
we return to Eq.(|67D we can see that the integral has a lower bound of R\/(f3 + 2) and this 
is achieved when (p)yim — 0- In this case we get a bound on 2m\jR\ given by 



, / 2mi 2mi 



(69) 



Thus we get the following bound 



'1 + 



09 + 2) 



(70) 



The approximation given in Eq. (|7"0|) holds only in the limit as f3 becomes large. In general, 
we can show that the expression in Eq. (|70|) is always less than 1 and monotonically increases 
with (3. For example, when j3 = we get 2m\jR\ < 0.944. 

We know that < (p)mn < (p)i- For any fixed value of /3, if (p)ma ~ (p)i the bound we 
get on 2ra\jR\ agrees with the monotonic bound as given implicitly by Eq.fl66|). As (p)ma 
reduces below (p) 1 the bound on 2m\jR\ increases monotonically and in the limit (p)Min ~ 
it reaches the bound given in Eq. ([7C|) and so is always strictly bounded away from 1. 



VII. CONCLUSIONS 

We have examined the ratio of the quasi- local mass to the circumferential radius, 2m/ R, 
for physically reasonable spherically symmetric isolated static configurations in general rel- 
ativity. We have demonstrated how the theory always places an upper bound on this ratio 
which lies strictly below the value it assumes when a horizon forms. This extends consider- 
ably earlier work on this question. 

The bounds we have derived do not take into account the stability of these static equi- 
libria. It would be interesting to know if and how these bounds get tightened when only 
stable configurations are considered. 

In || Morris and Thorne addressed the problem of constructing a static wormhole. If 
one has a spherical static wormhole one must have a minimal surface and thus p + Sc < 
somewhere. This raises an interesting question. Assume one has a spherical static wormhole 
and assume p > 0. How much 'exotic material' violating the strong energy condition does 
one need? 
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